Inflationary spectra and violations of Bell inequalities 
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In spite of the macroscopic character of the primordial fluctuations, the standard inflationary 
distribution (that obtained using linear mode equations) exhibits inherently quantum properties, 
that is, properties which cannot be mimicked by any stochastic distribution. This is demonstrated 
by a Gedanken experiment for which certain Bell inequalities are violated. These violations are in 
principle measurable because, unlike for Hawking radiation from black holes, in inflationary cosmol- 
ogy we can have access to both members of correlated pairs of modes delivered in the same state. 
We then compute the effect of decoherence and show that the violations persist provided the deco- 
herence level (and thus the entropy) lies below a certain non-vanishing threshold. Moreover, there 
exists a higher threshold above which no violation of any Bell inequality can occur. In this regime, 
the distributions are "separable" and can be interpreted as stochastic ensembles of fluctuations. 
Unfortunately, the precision which is required to have access to the quantum properties is so high 
that, in practice, an observational verification seems excluded. 



The inflationary paradigm successfully accounts for 
the properties of primordial spectra revealed by the com- 
bined analysis of CMBR temperature anisotropy and 
Large Scale Structure spectra |2|. In particular, it pre- 
dicts that the distribution of primordial fluctuations 
is homogeneous, isotropic and Gaussian, and that the 
power spectrum is nearly scale invariant (simply because 
the Hubble radius was slowly varying during inflation). 

Surprisingly, inflation implies that density fluctuations 
arise from the amplification of vacuum fluctuations Q; 
because of backreaction effects, the vacuum is indeed the 
only possible initial state Q . In addition of being ampli- 
fied, the modes of opposite wave-vectors k and — k end 
up highly correlated. More precisely, using linear mode 
equations, the vacuum evolves into a product of two- 
mode squeezed states [M IE • The highly squeezed 
character of the distribution implies the vanishing of the 
variance in one direction in phase space. This direction 
is that of the decaying mode ■ The observational con- 
sequence of this squeezing are the acoustic peaks in the 
temperature anisotropy spectrum 0, ll(Tj | . 

In spite of the macroscopic character of the mode am- 
plitudes, wc shall show that the inflationary distribution 
is still entangled in a quantum mechanical sense. To 
prove this, we shall provide observables able to distin- 
guish quantum correlations from stochastic correlations. 
At this point, it is important to notice that, unlike for 
Hawking radiation from black holes, we have in principle 
access to the purity of the state since, both members of 
two-mode sectors in the same state can be simultaneously 
observed on the last scattering surface 0] ■ 

Another important element should now be discussed: 
the linear mode equation is only approximate. Indeed, 
even in the simplest inflationary models there exists grav- 
itational interactions which couple sectors with differ- 
ent k's, and induce non-Gaussianities |l2j |. However, as 



in the BCS description of super-conductivity |lj|, the 
weakness of the interactions allows to approximate the 
distribution by a product of Gaussian two-mode distri- 
butions UJjhJ. The non-linearities will then affect the 
power spectrum as if some decoherence effectively oc- 
curred. In this sense, inflationary distributions belong 
to the class of Gaussian homogeneous distributions ob- 
tained by slightly decohering the standard distribution 
derived with linear mode equations. Notice also that in 
general, we have an experimental access to the state of a 
system only through a truncated hierarchy of it's Green 
functions, the Gaussian ansatz being the lowest order 
(Hartrcc) approximation. 

In the absence of a clear evaluation of the importance 
of non-linearities |25j . it is of value to phenomenologi- 
cally analyze the above class. It is characterized by three 
fc-dependent parameters. The first governs the power, 
see nk in The second gives the orientation of the 
squeezed direction in phase space, whereas the third con- 
trols the strength of the correlations between modes with 
opposite momenta. The latter is strongly affected by de- 
coherence effects, and shall be used to parameterize the 
decoherence level. It has been understood 0, ^3 that 
this level cannot be too high so as to preserve the well 
denned character of the acoustic peaks. However what 
is lacking in the literature concerning the quantum-to- 
classical transition is an operational identification of the 
subset of distributions exhibiting quantum correlations. 

To fill the gap, we propose a Gedanken experiment 
which shows that certain Bell inequalities are violated 
when using the standard distribution. We then show 
that the violation persists provided that the decoherence 
level lies below a certain threshold. Finally we point that 
there exists a higher threshold above which no violation 
of any Bell inequality can occur. The corresponding dis- 
tributions are separable (see below for the definition) and 
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can be interpreted as stochastic ensembles. 

In inflationary models based on one inflaton field, the 
linear metric (scalar and tensor) perturbations around 
the homogeneous background are governed by massless 
minimally coupled scalar fields obeying canonical com- 
mutation relations 01 ■ The scalar metric perturbations 
are driven by the inflaton fluctuations and correspond 
to perturbations along the background trajectory, called 
adiabatic perturbations ^3] • At the end of inflation, the 
homogeneous inflaton condensate decays and heats up 
matter fields. After inflation, during the radiation domi- 
nated era, the adiabatic perturbations correspond to den- 
sity perturbations of the matter fields (radiation, dark 
matter, ...) which all start to oscillate in phase. The fluc- 
tuations orthogonal to these, called iso-curvature, are not 
excited on cosmological scales in one inflaton field mod- 
els. Therefore, in the linear approximation, the phase 
and amplitude of the k-th Fourier mode of each mat- 
ter density fluctuation is related, via a time dependent 
transfer matrix, to the value of 0k and its time deriva- 
tive evaluated at the end of inflation (0 being the canoni- 
cal field governing scalar metric fluctuations during infla- 
tion) . This implies that the properties of the correlations 
of the density fluctuations are the same as those of 0k- 

We now briefly outline how one obtains highly squeezed 
two- mode states 0,0- During inflation, in the linearized 
treatment, each 0k evolves under its own Hamiltonian 
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where r\ is the conformal time dr\ = dt/a and a is the scale 
factor. To follow the mode evolution after the reheating 
time rj r , we continuously extend the inflationary law to 
a radiation dominated phase wherein a oc r\. In quantum 
settings, the initial state of the relevant modes (i.e. today 
observable in the CMBR) is fixed by the kinematics of 
inflation : these were in their ground state about 70 e- 
folds before the end of inflation (the minimal duration of 
inflation to include today's Hubble scale inside a causal 
patch). From horizon crossing k/a = H till the reheating 
time, (k 2 — d 2 a/a) in is negative. As a result, at the 
end of inflation, the initial vacuum has evolved into a 
tensor product of highly squeezed two-mode states. 

The resulting distribution belongs to the class of Gaus- 
sian homogeneous distributions, see [LH for more details. 
These are characterized by their two-point functions, best 
expressed as 



(4 a v ) =n k <5 3 (k-k'), 



(a k a k ,) = c k 5 3 (k + k'). (2) 
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The destruction operator a k is defined by a k e 
y/k/2(cj>k + idri^k/k) where 0k is evaluated at r] r . The 
mean occupation number governs the power spectrum, 
as shall be explained after Eq. (0. To meet the ob- 
served r.m.s. amplitude of the order of 10 -5 , one needs 
rik ~ 10 100 , i.e. highly excited states. The phase arg(cfe) 



gives the orientation of the squeezed direction in phase 
space at ry r . In inflation, using the above phase con- 
ventions, one gets arg(cfc) = n + 0(n k ). Finally, the 
norm of c k governs the strength of the correlations be- 
tween partner modes k, — k, i.e., the level of the coherence 
of the distribution. To parameterize the (de)coherence 
level, we shall work at fixed n and arg(c) (in the sequel 
we drop the k indexes), and write the norm \c\ as 



|c| 2 = (n + l)(n-£). 



(3) 



The standard distribution obtained in the linear treat- 
ment is maximally coherent and corresponds to 5 = 0. 
The least coherent distribution, a product of two ther- 
mal density matrices, corresponds to 5 = n. 

The physical meaning of 5 is revealed by decompos- 
ing the adiabatic modes in terms of the amplitudes (g, d) 
of the growing and decaying solutions. Taking into ac- 
count the time dependence of the corresponding transfer 
matrix, any matter density fluctuation can be used. For 
simplicity, we shall use the massless field extended in 
the radiation dominated era. In this case, the transfer 
matrix of ak is simply e~ %kr] . Decomposing 



4>k(v) = 5k 



sin(fc?7) 5 cos(krj) 



Vk 



Vk 



(4) 



Eqs. J2J gi vc 



(5k5 k > = « + \ - Mc) = 2n (l + 0(^j , 
(d k dl) = n+ l - + Re(c)= 5 - + 0{n-^). 



(5) 



The last expression in each line is valid when the deco- 
herence is weak, i.e. i5 <C n. In this regime, the power 
spectrum P k = fc 3 (0k(??)0-k(??)) — k 2 rik s\n 2 (krj) is dom- 
inated by the growing mode. At fixed rj, it therefore dis- 
plays peaks and zeros as k varies. From the last equation 
|[SJ , one sees that the decoherence level S fixes the power 
of the decaying mode. (The same conclusions would have 
been reached had we considered dark matter or temper- 
ature perturbations.) 

Even though Eqs. @ univocally determine the corre- 
sponding (Gaussian) distribution, they are unable to sort 
out the distributions possessing quantum properties from 
those which have lost them, or in other words, to deter- 
mine the ranges of <5 characterizing these two classes. To 
operationally do so, it is necessary to introduce operators 
which are not polynomial in g k and d k |2(| . 

In what follows, we shall use operators based on 
coherent states. These obey Ok|«,k) = v\v, k) and 
a_k|u>, — k) = w\w, — k). They are minimal uncertainty 
states and each of them can be considered as the quan- 
tum counterpart of a point in phase space, here a classi- 
cal fluctuation with definite phase and amplitude. This 
correspondence is excellent in the regime n» 1. More- 
over, they play a key role when considering decoherence: 
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when modes are weakly coupled to an environment, the 
reduced density matrix becomes diagonal in the basis of 
coherent states or other minimal uncertainty states 

Coherent states are particularly useful in our context 
because they will allow us to sort out entangled quan- 
tum distributions from stochastic ones. The reason is 
that coherent states can probe the detailed properties of 
the distribution. In particular, the probability to find a 
particular classical fluctuation is given by the expectation 
value of the projector on the corresponding (two-mode) 
coherent state, namely 



H(v, w) = \v, k) (v, k| <g> \w, — k) (w, — k| 
The probability is 

Q(v,w;S) = Tr[p 2 (6)IL(v,w)], 



(0) 
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n+1 
1 



exp 



(n + I) 



l + <5 



exp 



\w — w(v)\* 
l + S 



(7) 



where pviS) is the matrix density of the two- mode sys- 
tem. We have written Q(v,w;S) in an asymmetric form 
to make explicit the power of the growing mode (= n+1), 
and the much smaller width (=1 + 5) governing the dis- 
persion of the values of w around w(v) = v*c/(n + 1), 
the conditional amplitude of the partner mode, given v. 
Had we used a projector on a one-mode coherent state, 
we would have gotten only the first Gaussian. In fact, 
as we shall see, to have access to the (residual) quan- 
tum properties of the distribution, one must use the two- 
mode projectors @). As explained in ^lj, these pro- 
jectors also allow to compute conditional values which 
cannot be expressed in terms of mean values. For in- 
stance, Tr[ITpn <fi(n, x)] gives the space-time pattern of 
fluctuations when the set of configurations specified by 
the projector II is realized. 

Given the macroscopic character of mode amplitudes 
in inflationary cosmology, it is remarkable that the pro- 
jectors (JBJ can violate Bell inequalities. To understand 
the origin of this possibility, it is necessary to define the 
class of separable states 20] . A two-mode state is said 
separable if it can be written as a positive sum of prod- 
ucts of one-mode density matrices Separable Gaussian 
states can all be written in terms of the projectors as 



P2 p -(S)= f— — P(v,w;S)IL{v,w)- (8) 



,(9) 
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The function P is given by 



P{v, w; 5) 



1 



-exp 
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x exp 


A'/n 



with w = cv*/n and A' — n 2 — \c\ > 0. The latter 
implies |c| < or 8 > n/(n + 1) ~ 1 for n 3> 1. (The 
limiting case \c\ = n, 8 = n/(n + 1) is interesting: the 
second exponential becomes a double Dirac delta which 
enforces w = cv*/n = ~v* in phase and amplitude. In 
other words, for each two-mode sector, there is only one 
fluctuating quantity, since the second mode is completely 
fixed by its partner. In inflationary cosmology, the corre- 
sponding density matrix can be viewed as the quantum 
analogue of the usual stochastic distribution of growing 
modes. Indeed, the entropy of this quantum distribu- 
tion is ln(n) per two-mode, and this is the entropy of the 
stochastic distribution for each growing mode . This 
quantum-to-classical correspondence is corroborated by 
the fact that off-diagonal matrix elements of p(S) in the 
coherent state basis vanish precisely when S > 1. 

The physical meaning of separable states comes from 
the fact that all states of the form iJHJl can be obtained by 
the following classical protocol ji^l : when a random gen- 
erator produces the four real numbers encoded in (v,w) 
with probability P, two space-like separated observers 
performing separate measurements on the subsystems k 
and — k respectively, prepare them into the two-mode 
coherent state Non-separable states can only be 

produced by letting the two parts of the system interact. 
Only these are quantum mechanically entangled. 

By construction, the statistical properties of separable 
states can be interpreted classically. In particular, they 
cannot violate Bell inequalities [20j. In what follows we 

" be- 



(10) 



shall study the "Clauser-Horne" inequality |2JJ, 122 
cause it is based on Q of Q. It reads 

C{v,w;5) = [Q{0,0;S) + Q{v,0;6) 



+ Q(0, w; S) - Q(v, w; 5)} x ( ] 



We can now search for distributions, i.e. values of 5, 
and for configurations v and w which maximize C. The 
maximization with respect to w gives axg(c*vw) = tt and 
\w\ = \v\. We fix the arbitrary phase of v by 2arg(w) = 
arg(c), so that C is maximum along the 'line' w = —v. In 
Fig. 1 we have plotted C(v, —v, 5) for three values of S. 
The maximum with respect to the norm of v is reached 
for 



\vM(S)f 
l + S 



In 


1 + 


/n-S 









1 + 2, 



In 2 

1T 



1 + 0(5/")] -(H) 



The maximal value is 
1 



Cm (8) 



1 + 2^ + °hr 



2{l + 8) 

The inequality (|10f) is thus violated for 



5<- 1+ y 24/3 ^ 0.095, 



(12) 



(13) 
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FIG. 1: The loss of violation as decoherence increases. We 
have represented C(v, —v, 5) as a function of x = |ti| 2 for n = 
100 and for three values of 5: (upper), 0.05 (middle), and 
0.1 (lower). The horizontal line (C = 1) is the maximal value 
allowed by classically correlated states. 

irrespectively of the value of n when n» 1. 

From the last two equations we learn that Bell inequal- 
ity \1U\) is violated by the standard inflationary distribu- 
tion (5 — 0). Notice that this violation is maximal, as 
one might have expected, since the two-mode correla- 
tions are the strongest in this state. More importantly, if 
S obeys (|X3I) , the violation persists in the regime of highly 
amplified modes obtained in inflationary cosmology. 

In conclusion, the principle results of this Rapid Com- 
munication are the following. First, in spite of the 
macroscopic character of adiabatic fluctuations, the stan- 
dard inflationary distribution possesses quantum features 
which cannot be mimicked by any stochastic distribu- 
tion. Second, these features arc operationally revealed 
by a well defined procedure based on the violation of the 
Bell inequality I|1U|) . Third, the projectors used in this in- 
equality have a clear meaning in cosmology: they give the 
probability that a particular semi-classical fluctuation be 
realized. Fourth, the mere existence of decoherence ef- 
fects is not sufficient to eliminate the quantum properties. 
To do so, decoherence should be strong enough so as to 
induce S > 1, that is, so that the distribution becomes 
separable. 

The threshold value (5=1 therefore plays a double 
role. First, as previously noticed, the distribution with 
5 = 1 possesses an entropy (= Inn per two-mode) which 
is equal to that of the classical distribution of growing 
modes. Second, separability is the condition for distin- 
guishing quantum from classical distributions, see e.g. 
[23j where it was used to define the time of decoherence. 
To our knowledge, besides the present work, this crite- 
rion of the study of the quantum-to-classical transition 
has not been used in inflationary cosmology. 

Let us now briefly address two additional questions. 
Firstly, to what extend the violation of the inequality 
(|10|l is verifiable ? We start by pointing out that there is 
no physical principle which prevents evaluating the four 
terms in Eq. (|1(J|> . Because of isotropy, in a given co- 
moving volume (e.g. a sphere of radius R), we have, 



for a given wave vector norm k = |k|, about (kR) 2 adia- 
batic modes all characterized by the same two-mode den- 
sity matrix. This is true before and after the reheating, 
and also irrespectively of the decoherence level. Finally 
this is still true when considering the projection of the 
adiabatic modes on the last scattering surface. Indeed, 
for sufficiently high angular momentum, there exist an 
ensemble of well aligned two-modes with both members 
living on the last scattering surface O ne can thus 

accumulate statistics to measure the four observables of 
Eq. I|1U|) . Unfortunately, an observational verification of 
the inequality Eq. (|10|) seems excluded. Indeed the cos- 
mic variance, which is of the order of the mean amplitude 
(hence proportional to n ~ 10 100 ), is much larger than 
the required precision, which is given by the spread of 
the coherent states (=1) H3 

The second question concerns the value of 5 in realistic 
inflationary models. This interesting question deserves 
further study, see the Added Note. Let us here sim- 
ply compare the critical value 6 = 1 separating quantum 
from stochastic distributions to the expected level of non- 
Gaussianitics. At the end of inflation, the two-point func- 
tion of the gravitational potential ^ is conventionally |3| 
parameterized by the coefficient f^L entering the field 
redefinition * = [cf> + f NL (f> 2 /(aM )] /(y/eMa) where <f> 
is our Gaussian field during inflation, M is Planck mass, 
and e the slow roll parameter. It has been observationally 
limited to —58 < Jnl < 134 Q, while theoretical calcu- 
lations give fNL = O(10~ 2 ) for the inflationary phase. 
On one hand, the variation of the power spectrum of "J 
is therefore AP/P ~ /jvL-P where P is the power spec- 
trum in the linear approximation (~ 10~ 10 ). On the 
other hand, using JSJ, one gets AP/P = S/n. Therefore 
f NL = 1(T 2 corresponds t S ~ nPf% L ~ 10 86 . This 
indicates that the minimal source of decoherence, the 
non-linear interactions during inflation, should be strong 
enough to give rise to separable distributions. 
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